Some exact results for Boltzmann's annihilation dynamics 



Frangois Coppex/ Michel Droz,^ Jaroslaw Piasecki,^ Emmanuel Trizac,"^ and Peter Wittwer-'^ 

^Department of Physics, University of Geneve, CH-1211 Geneve 4, Switzerland 
. . . ^Institute of Theoretical Physics, University of Warsaw, PL-00 681 Warsaw, Poland 

' Laboratoire de Physique Theorique, Batiment 210, Universite de Paris-Sud, 91405 Orsay, France 

o ; 

, The problem of ballistic annihilation for a spatially homogeneous system is revisited within Boltz- 

. mann's kinetic theory in two and three dimensions. Analytical results are derived for the time 

evolution of the particle density for some isotropic discrete bimodal velocity modulus distributions. 
According to the allowed values of the velocity modulus, different behaviors are obtained: power law 
decay with non-universal exponents depending continuously upon the ratio of the two velocities, or 
exponential decay. When one of the two velocities is equal to zero, the model describes the problem 
' of ballistic annihilation in presence of static traps. The analytical predictions are shown to be in 

agreement with the results of two-dimensional molecular dynamics simulations. 
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Q . I. INTRODUCTION 

In ballistically controlled reactions, particles with a given initial velocity distribution move freely (ballistic motion) 

^ in a d-dimensional space. When two of them meet, they annihilate and disappear from the system. This apparently 

Ci simple problem has attracted a lot of attention during the past years 

El II II Si, Ell 11 ills in III mill for 

c/2 the following reasons. First, this is one of the few problems of nonequilibrium statistical physics that can be exactly 
^ , solved in some cases, and second it models some growth and coarsening processes |15( . 

^ ■ This field was entered with the pioneering work by Elskens and Frisch 0, where a one-dimensional system with 
^ , only two possible velocities +c or — c was studied. Using combinatorial analysis, they showed that the density of 
'. ■ particles was decreasing according to a power law in the case of a symmetric initial velocity distribution. 

The investigation of this one-dimensional problem was generalized by Droz et al. to the three- velocity case where 
Q i the initial velocity distribution is given by tf(v]Q) = P+5{v — c) + pq5{v) + p^5{v + c) with p+ = p_ (symmetric 
Q ' case) and p+ -I- po + = 1 • It turns out that the decay of the particle density depends on the details of the initial 
velocity distribution. The following analytical results were obtained. For po < 1/4, the density n{v]t) of particles 
' with velocity v — {0,+c, — c}, behaves in the long-time limit as n{0;t) ~ t~^, n{±c;t) ^ t^^^^. When pq = 1/4, 
] n{0;t) ^ n{±c;t) ^ t'^^^. Finally, for po > 1/4, one finds that n{0;t) saturates to a nonzero stationary value, while 
\^ ' n(±c; t) decays faster than a power law. Moreover, it was shown that in one dimension, annihilation dynamics creates 
[ strong correlations between the velocities of colliding particles, which excludes a Boltzmann-like approximation. Pairs 
CO . of nearest neighbor particles have the tendency to align their velocities and propagate in the same direction |^. 

' An analytical investigation of the one-dimensional case with a continuous velocity distribution is much more difHcult. 
A dynamical scaling theory, whose validity was supported by extensive numerical simulations for several velocity 
distributions, led Rey et al. 7\ to the conjecture that all the continuous velocity distributions ^(v) that are symmetric, 
regular and, such that Lp{0) ^ are attracted in the long-time regime towards the same Gaussian-like distribution 
■ and thus belong to the same universality class. 

, For higher dimensions, most of the studies are based on an uncontrolled Boltzmann-like description H H or 
I ■ numerical simulations [TJ. However, based on phenomenological mean-field-like arguments, Krapivsky el al. have 
' studied the annihilation kinematics of a bimodal velocity modulus distribution in d > 2 dimensions In the case 
^ , of a mixture of moving and motionless particles they showed that the stationary particles always persist, while the 
density of moving particles decays exponentially. This approach contains unknown phenomenological parameters, and 
thus a complete comparison with the results obtained by numerical simulation is not possible. 

In a recent paper, Piasecki et al. |l4j gave an analytical derivation of the hierarchy equations obeyed by the reduced 
distributions for the annihilation dynamics. In dimension d > 1 for a spatially homogeneous system, and in the 
limit (the so-called Grad limit) for which the particle diameter a —^ 0, and the particle density n{t) — > oo such that 
I n{t)a'^^^ = A^^, where A is the mean free path, the hierarchy reduces to the Boltzmann-like hierarchy. This hierarchy 
propagates the factorization of the reduced fc-particle distribution in terms of one-particle distribution functions. Thus, 
if the initial state is factorized, the whole hierarchy reduces to one nonlinear equation for the one-particle Boltzmann 
distribution. For annihilation kinetics, the ratio of particle diameter to mean free path vanishes in the long-time limit 
and the situation becomes similar to the Grad limit discussed above for A — > oo. Thus the long-time limit of the 
annihilation dynamics (for d > 1) is likely to be adequately described by the nonlinear Boltzmann equation. 

A scaling analysis of the nonlinear Boltzmann equation led to analytical expressions for the exponents describin g th e 
decay of the particle density and of the root-mean-square velocity in the case of continuous velocity distributions 14] . 
In view of the different behaviors observed in one dimension for discrete or continuous velocity distributions, it is 
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relevant to study the case of distributions with discrete modulus spectrum in dimensions higher than 1. The goal 
of this paper is to investigate simple examples of this kind in three dimensions for which the non linear Boltzmann 
equation derived in Ref. jlj| can be exactly solved. The generalization of this approach to an arbitrary dimension is 
straightforward, and for the sake of comparison with numerical simulations, we shall also consider the two-dimensional 
situation in some detail. 

The validity of the Boltzmann description in the long-time limit will be confirmed by comparing our analytical 
predictions with the results obtained by a molecular dynamics simulation. 

The paper is organized as follows. In Sect. ^ we define the model. In Sect. IIIll the three-dimensional Boltzmann 
equation is solved analytically for a two velocity modulus (ci and C2) isotropic distribution. For simplicity we first 
consider the one velocity model Ci = C2 > that allows to draw interesting comparisons with the same model in one 
dimension. Then the implicit solution for the particle densities in the general case Ci > C2 > is established. It is 
shown analytically that in the long-time limit the particle densities decay according to power laws, with exponents 
depending continuously on the value of the velocity modulus ratio. We also find upper and lower bounds to the 
particle densities that are compared with the numerical solution of the dynamical equation. The particular case 
of a mixture of moving (c2 > 0) and motionless (ci = 0) particles is also investigated. It turns out that the 
particle densities decay exponentially to zero for the moving particles, and to a nonzero value for the motionless ones. 
This phenomenology is independent of space dimension, and in Sect. IIVI it will be shown explicitly to hold in two 
dimensions by implementing molecular dynamics simulations. This numerical method has the advantage of being free 
of the approximations underlying Boltzmann's dynamics and, therefore, provides an interesting test for the analytical 
predictions. Section Ivl contains our interpretations and conclusions. 



II. THE MODEL 



We consider a system made of spheres of diameter tr moving ballistically in three-dimensional space. If two particles 
touch each other, they annihilate and thus disappear from the system. We consider only two-body collisions. The 
initial spatial distribution of particles is supposed to be uniform, therefore it remains uniform during the evolution. 
Indeed, if the state is initially translationally invariant, then the free evolution preserves this symmetry. Annihilation 
dynamics adds the effect of binary collisions that depends only on the distance between particles, thus preserving 
the translational invariance. Finally, existing numerical simulations seem to be compatible with this assumption of 
homogeneity [T^ . We are interested in the time evolution of the number density of particles with a given velocity 
modulus. 

Let /i(v; t) be the distribution function of the density of particles in R'^ with velocity v e R'^ at time t. For spatially 
homogeneous states, the distribution function has the form 

h{^;t)=n{t)ip{v;t), (1) 

where (^(v; t) is the velocity probability density. In the long-time limit, Piasecki et al. |l4j have shown that the 
hierarchy satisfied by the reduced distributions approached the Boltzmann hierarchy. If the initial state is factorized, 
the nonlinear Boltzmann equation provides then the complete description of annihilation dynamics. 



d_ 

dt 



fi{vi;t) = a^ [ d9 e{-9-vi2)i^-vi2) [ dv2 |vi2|/i(vi;i)/i(v2;t). (2) 
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Here 9 is the Heaviside function, V12 = vi — V2 the relative velocity of two particles, V12 — V12/U12 a unit vector, 
V12 = |vi2|, and the integration with respect to da is the angular integration over the solid angle. 

We consider spherically symmetric initial conditions /i(v;0), v — |v|. This symmetry property is propagated by 
the dynamics. The Boltzmann equation jSJ then takes the form 



{u + vY — |u — v\^ 



(3) 



Equation Q is a nonlinear homogeneous integral equation for the distribution function t). A simplification arises 
if the initial velocity distribution has a discrete modulus spectrum. This spectrum is preserved by the annihilation 
dynamics as no new velocities are created. A simple case is provided by the bimodal distribution 

^(«' 0) = -^Kv - ci) + ^^S{v - C2), (4) 

where C2 > ci > and A denotes the fraction of particles with velocity modulus ci. 



3 



III. EXACT RESULTS 

Before addressing the general case, we first consider the single-species problem where C2 — ci > 0. 

A. Single-velocity modulus distribution 

Setting C2 — ci — c > in Eq. one obtains from Eq. 



From the kinetic equation 0, we find 



whose solution is 



fi{v;t)^n{t)^S{v~c). (5) 



|n(^) = -|w2cn(^)^ (6) 



f + ^na'^noct 



where nlO) ~ uq. A striking observation is that, in the limit t 00, the density (0 becomes independent of its initial 
value no- Note that the same phenomenon is also present for simple diffusion limited annihilation such as A + A ^ 0, 
when the dimension of the system is larger than 2 . 

Contrary to the one-dimensional case for which it has been rigorously shown that the density decays proportionally 
to one sees from Eq. (0 that in three dimensions, Boltzmann's dynamics is faster as the density decays 

according to t~^, which is the mean-field value Q. We note, however, that the same behavior n(t) oc 1/t holds 
in all dimensions within Boltzmann's kinetic theory (and in fact, more generally within the framework of a scaling 
analysis of the hierarchy governing the dynamics of ballistic annihilation 14]). This discrepancy between Boltmann's 
prediction and the exact result in one dimension illustrates the crucial importance of dynamical correlations when 
d = 1. On the other hand, as suggested in Ref. |14| and explicitly shown below by molecular dynamics simulations, 
the nonlinear Boltzmann equation is relevant for describing the long-time dynamics of ballistic annihilation when 
d > 2. In this case the particles are very diluted and no dynamical correlations can develop during the time evolution, 
which would violate the molecular chaos hypothesis. 

B. Mixture of particles with two nonzero velocity moduli 

Consider the case where particles with velocities Ci > and C2 > Ci are initially present. Thus fi{v;t) is of the 
form 

Mv;t) = X{t)^^S{v - ci) + Y{t)^^6{v - C2), (8) 

where X{t) and Y{t) are, respectively, the densities of particles with velocities ci and C2. They add up to the total 
density X{t) + Y{t) — n{t). Upon rescaling the time according to r = te27r(T^/3, it follows from Eq. Q that 

X{t) = -4jX{Tf -{3 + j^)X{t)Y{t), (9a) 
y(r) -4y(r)2 - (3 + 7^)X(r)y(r), (9b) 

where < 7 = Ci/c2 < 1, and the overdot denotes time derivative with respect to r. 

The set of equations © is a nonlinear homogeneous system of coupled differential equations with constant coeffi- 
cients. An implicit solution can be obtained by introducing the function V{t) defined as V{t) — Y{t)/X{t). From 
Eq. (O we get 

dY ^ W^ + i3 + -f^)V 

dX 47 + (3 + 72)T/ ' ^ ^ 

so that 

dX 47 (3-^7^)1/ 



X (1 - 72)1/2 (3 



dV. (11) 
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Integrating Eq. yields 




(12) 



with a = 47/[(l - 7)(3 - 7)] > 0, /3 = (3 + 7')/(l - 7') - " > 0, ^(0) - = l^o/^o, ^0 = ^(0), Yo = Y{Q). 
The special case of 7 = will be discussed in Sec. IIII CI hence from now on we assume that 7 > 0, so that a > 0. 
Equations (O can also be written as 

Multiplying the right-hand side (RHS) of Eq. I|13af) by Xq and equating it with the derivative of the RHS of Eq. l(T^ . 
one obtains upon integration from to r the relation 



Xqt = / du 



^« lid 



y 47 + (3 + 72)u I du 




(14) 



Equation p4(l implicitly defines the time dependence of the function V{t). The procedure to obtain the densities X{t) 
and Y{t) from Eq. (|14|l is as follows. The integration in Eq. (|14|) leads to Appel functions, that may be inverted (at 
least numerically) in order to give V{t). The insertion of V{t) in Eq. (|12l) then gives X{t). It is then straightforward 
to obtain Y{t), having determined V{t) and X{t). The structure of the implicit relation H14(l permits us to establish 
interesting analytical results. 

First, let us investigate the long time behavior of the particle densities X{t) and Y{t). When t — > 00, the LHS of 
Eq. H14II diverges linearly which implies that limT^oo Vi^) — 0. So, in the long-time limit, the implicit relation 114|) 
leads to the asymptotic formula 



^or'.-i(^) + , (15) 



A-i \V J V3-7 
so that 

y(r) ™ K, (^^^0 + 1 j (47Xor)-i/". (16) 

As lim^^o y{T) = 0, Eq. (|13a|l takes the asymptotic form d/dT{\/ X) = Aj. Hence we conclude that 

X{r) 3^,-1. (17) 
On the other hand, from Eqs. Ijl3b|) and (|16|l we find the long-time relation 

' ^ ' ^ (3 + 7^)^ f l^Vo + l) (47Xor)^/" , (18) 



dr \Y J ' ' 'Vb V3-7 

which upon integration yields 

r(r) ^(47X0)-^" (}±2v, + iy^%-(3+7^)/47. (19) 
47 V3-7 / 

Note that the exponent for the density Y{t) is a function of the ratio 7 = ci/c2 and thus is nonuniversal. In the limit 
7^1 one recovers the asymptotic behavior of the single- velocity modulus distribution (see Sect. IIII^ . 

Second, we may find analytical upper and lower bounds for X(t) and Y{t). Granted that a > and /3 > 0, 
the integrand of Eq. H14|l is a strictly monotonic decreasing positive function of u, therefore V{t) < Vq for all 
T > 0. Considering that (47)^^ > [47 4- (3 + 7^)11]^^ for u > 0, the insertion of Eq. (fT^ in Eq. ((TI|l provides the 



5 



inequality Xqt < {Xq/X — 1) /4j, which leads to an upper bound for X{t). On the other hand, the inequality 
(47 + (3 + j'^)u)^^ < [47 + (3 + 7^)Vo]^"'^ yields a lower bound, so that we finally get 



^0 



1 + [47X0 + (3 + j^)Yo]t 



< Xir) < 



1 + 47X0T 



Note that for times such that 



47Xor > 1, 



(20) 



(21) 



the upper bound 12011 coincides with the exact asymptotic relation (|17|l . The same kind of analysis as that leading to 
Eq. Iini) yields the upper bound. 



< Y{t) < 



Yn 



1 + (3 + 72)Xor 



(22) 



The width defined by the difference of the bounds in both cases (|20|l and (|22|l is O (t^^). Figures ^ and |21 show the 
numerical solution for X{t), Y{t), the bounds (|20|l and (|22|l . as well as their asymptotic behaviors H17|) and IjlOf) on 
a logarithmic scale. 
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FIG. 1; Upper and lower bounds 12011 (dotted lines) as well as the numerical solution of the set of equations @ for X{t) with 
Xo — Yq, 7 = 0.2 (continuous line). The inner logarithmic plot shows indeed the power law behavior X{t) ~ for r 00, 
where the asymptotic solution H17|l is represented by the dashed straight line. Moreover, in this regime the solution converges 
to the upper bound iCT . 



The knowledge of the numerical solution (see Figs. ^ and [2J) allows to determine the crossover time, separating the 
early and long-time (power law) regimes. 



C. Mixture of moving and motionless particles 

We now consider a particular case of Sec. 1111 lH that we solve exactly in the asymptotic limit r — > 00. The system is 
now characterized by a certain number of motionless particles (zero velocity, ci = 0) whereas the rest of the particles 
have a given nonzero velocity modulus. Thus, setting 7 = in Eq. (|12|l . inverting the relation in order to find 
V — V{X), then making use of Eq. (|13a|l with 7 = leads to 

d f l\ Yo \ f X^'^^^ 
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FIG. 2: Upper bound p2|l (dashed lines) as well as the numerical solution of the set of equations Q for Y{t) with Xq = Yq, 
7 = 0.2 (continuous line). The inner logarithmic plot of the numerical solution shows indeed the power law behavior Y{t) ~ 
^-(3+7 )/47 £qj- j- where the asymptotic solution 11911 is represented by the dashed straight. Furthermore, the use of both 

the upper bound 11711 and the asymptotic form 1191 allows to find an analytical approximation for Y{t), which turns out to be 
exact in both limits t and t — > oo. 



The integration of Eq. 123|) yields 
3XoT - 1 



3[i-^] + f i-{j^y 



In 



a-l 



a-(Xo/X)i/3 



(24) 



with a = 1 + Yq/'SXq. In the asymptotic hmit t ^ oo, the LHS of Eq. H24|l tends to +oo. The density X{t) cannot 
thus tend to zero, and must approach a strictly positive value X{oo) = X^c > 0. For t —>■ oo, all terms on the RHS 
of Eq. (|24|l but the logarithmic one approach a finite limit. This implies the asymptotic behavior 



X^ = lim X{t) 



Xn 



Xo 



3Xn 



> 0. 



(25) 



a- \3Xo + Yo, 

Replacing X by Xoo in all terms of Eq. 12411 except the logarithmic one and then inverting the relation t{X), we find 

T — i- OO Q T- — > OO 

X{r) ~ Xoo[l-e(Xo,ro;r)]"' ~ + 3e(Xo, Fo; r)], (26) 



where 



and n = 1/3 + a/2 



Hence we have 



.(^o,ro;r)^l^Lf^exp(-3X.r), 
lla^/16 < 0. Making use of Eq. ^ with 7 = 0, we find 



Yir) 



3X^e{Xo,Yo;T). 



X(t) 



X^ 



Y{t). 



(27) 



(28) 



(29) 



There is a qualitative difference from the case ci > 0. As shown in Fig.O the density of particles at rest approaches 
the asymptotic value Xoo > exponentially fast, while the density of moving particles goes to zero exponentially. 
Table U summarizes the long-time behavior for the different cases. 

Note that generalizing our results to any dimension d > 2 is straightforward (see below for the case d = 2). The 
algebraic or exponential decay of the particle densities hold irrespective of d. In particular, for the general case ci > 
the exponent of the density of "slow" particles is independent of d so that X{t) t^^. Finally the relation H29|) 
still holds. 
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12 3 4 

FIG. 3: Linear-logarithmic plot of the numerical solution of the set of equations ^ for Xq = Vo, 7 = (continuous lines). 
The asymptotic relations 12611 and 1281 are shown by the dashed lines, and the asymptotic limit (1251 by the dotted line. 



TABLE I: Summary of the long-time density behavior in three dimensions. 





C2 — Ci > 


C2 > Ci / 


C2 > Ci = 


X{r) 






Ji:oo[l + 3exp(-3Xoor)] 


Y{r) 




^-(3+7^)/47 


Xoo3exp(-3XooT) 



IV. COMPARISON WITH NUMERICAL SIMULATIONS 

The analytical predictions obtained in the preceding section rely on the validity of the molecular chaos assumption, 
leading to the Boltzmann equation. It is therefore instructive to compare these predictions to the results of molecular 
dynamics (MD) simulations, where the exact equations of motion of the particles are integrated (see Ref. jl^ for more 
details concerning the method). 

A. Analytical results in two dimensions 

MD simulations are most efficiently performed in two dimensions, where the best statistical accuracy can be 
achieved. We consequently repeat the analysis of Sees. ^ and IIIII for a two-dimensional system. Introducing the 
rescaled time r = 2'iTCfC2t, one obtains the counterpart of Eqs. in the form 

X{t) = -47X(t)2 - K(7)X(r)y(T), (30a) 
Y{t) = -AY(Tf -k{^)X{t)Y{t), (30b) 

where ^(7) — d(py/ 1 + 7^ — 27 cos ip. In the limit t ^ 00, the solution of the system H30|) reads 

Xir) ™ gi^r^i, (31a) 
r(r) ^i^lXor^^ (ifl^Vo + lY'r-/'^ (31b) 
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where n = i'y / {k — Aj) and v = k/ {A — n) — ji. On the other hand, taking the hmit 7 in Eqs. (|30|l and solving 
the corresponding system leads to the long-time behavior 



X[t) ^ Xoo[l+£2(^o,>^o;r)] 
Y{t) X^S2{Xo,Yo;t), 



(32a) 
(32b) 



where 



and X = 4/7r — 1. In Eqs. (|32|l . we have 



^0 (1 + Voxy/^ 



(33) 



£2(^0, Yo; t) = 7rV^o(^+^"^)"''""' exp {-JX^/Xo) exp {-^X^t) 



(34) 



and 



J = du \n{u) 
Jo 



(35) 



B. Molecular dynamics simulations 

MD simulations have been implemented with systems of typically N = 10^ to 4 x 10^ spheres in two dimensions 
(discs). Periodic boundary conditions were enforced, and low densities considered, in order to minimize the excluded 
volume effects discarded at the Boltzmann level (note that these effects are necessarily transient since the density 
decreases with time). 

Figure 21 compares the MD results obtained with 7 = 1/10 to the predictions of Eqs. H31|) (for 7 = 1/10, the time 
decay of the "fast" particles is governed by the exponent k/47 ~ 7.9). Although the large-time behaviors for X and 
Y are compatible with those given by Eqs. it may be observed that the corresponding asymptotic regime is 

difhcult to probe, even for large systems. The parametric plot [or "trajectory" 1^(X)] shown in the inset is however 
in agreement with the relation Y cx X"/*^^ deduced from Eq. %M\ . 

We have also performed MD simulations for a mixture of moving and motionless particles (7 — 0), where it is 
expected that the density X of particles at rest decreases down to a nonvanishing value Xoc- In the situation of an 
equimolar mixture {Xq = Yb)j we have Vb = 1 so that according to Eq. (|33|l . Xoo/Xq ~ 0.414. The MD simulations 
are in agreement with this scenario, and we find Xac/ X^ 2± 0.408 irrespective of the initial conditions for a system 
with initially N ~ 2 x 10^ particles. The results for the time dependence of X and Y are displayed in Fig. |S1 We 
conclude that the numerical simulations are again in agreement with the prediction of Boltzmann's kinetic theory. 



V. CONCLUSIONS 



We have shown that for some simple spatially homogeneous systems, characterized by a velocity distribution with 
a discrete velocity modulus spectrum, it is possible to find the exact solution for the nonlinear integral equation 
describing the dynamics of ballistic annihilation. These results, obtained at the level of a Boltzmann equation, have 
been validated by explicit comparison with molecular dynamics simulations in two dimensions. 

For a single- velocity modulus distribution, the particle density of the model decays asymptotically as n{t) t~^, 
irrespective of space dimension. It was however rigorously shown that in one dimension, the decay is slower, n(t) ^ 
^-i/2_ ^j^jg difference is a consequence of the fact that in one dimension strong dynamical correlations are created 0, 
which invalidate the approximation underlying Boltzmann's dynamics. In higher dimensions, the Boltzmann equation 
becomes exact in the long-time limit. 

In the case of a distribution with two different finite nonzero velocity moduli, we found that both particle densities 
decay for a large time according to a power law. The interesting feature is that the density of the slow particles decays 
as t~^, while the density of the fast particles decays more rapidly (e.g., as ^^(-^+7 )/47 \^ three dimensions), with a 
nonuniversal exponent depending continuously on the velocity modulus ratio 7 — cxjci. A rough criterion for the 
crossover time separating the short- and long-time regimes has been given in Eq. %2,\\ . 

Finally, the case ci = leads to a particularly interesting behavior. Independently of the initial conditions, the 
densities of the moving and the motionless particles both decrease exponentially fast; however down to a nonzero 
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FIG. 4: Log-log plot of the densities X (upper curve) and Y (lower curve) as a function of rescaled time, as obtained in the MD 
simulations of a two-dimensional system with 7 = 0.1. The initial condition corresponds to an equimolar mixture [Xq — Yo) of 
N — 2 X 10^ particles, with reduced density {Xo + Yo)a'^ — 0.1 at r = (both species have the same diameter). The dashed 
lines have slopes —1 and —7.9 [as predicted by Eqs. 1311 1. Inset: log-log plot of y as a function of X, where the broken line 
has slope —7.9. 

value for particles at rest. This behavior is quite different from that observed in the one-dimensional case where the 
initial value of the density of motionless particles plays an important role in the long-time regime. This difference 
between one dimension and higher dimensions reflects once again the important role played by the dynamically created 
correlations for d — 1. 

The case with motionless particles can be viewed as a problem of ballistic annihilation of particles with one- velocity 
moduli moving in a random medium containing immobile traps (the motionless particles) that can capture a moving 
particle and then disappear. Here again, the situation can be compared to similar problems in diffusion limited 
annihilation where the presence of traps can modify the long-time dynamics from a power law to an exponential 
decay [iTj . 

It would be interesting to compare the above theoretical predictions with some experimental data. Besides growth 
and coarsening problems, ballistic annihilation could model other physical systems such as, for example, the flu- 
orescence of laser excited gas atoms with quenching on contact [l^l- However, the correspondence between such 
experimental situations and our model is not yet close enough to allow comparison. We would be highly interested in 
the knowledge of other physical systems that could be described by the models studied here. 
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